AD-A208  111 


•  »» 


t 


NAVAL  POSTGRADUATE  SCHOOL 


Monterey,  California 


% 


* 


DTIC, 

FIECTE  I 
MAY  2  6  19891 


4 


UNCLASSIFIED 


SELURI T  Y  CLASSIFICATION  OF  TH'S  PAUL 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704  0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


a  SECURITY  CLASSIFICATION  AUTHOR'TY 


2b  DEClASSiFiCATION  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGAN'ZATlON  REPORT  NUMBER(S) 


1b  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUMBER'S) 


6a  NAME  OF  PERFORMING  ORGANIZATION  I  6b  OFFICE  SYMBOL  I  7a  NAME  OF  MONITORING  ORGANIZATION 

|  (if  applicable)  I 


Naval  Postcraduate  School!  Code  35 


6c  ADDRESS  f C/f>,  Srafe,  and  ZIP  Code) 

Monterey,  California  93943-5000 


Naval  Postqraduate  School 


7b  ADDRESS(Cify,  State,  and  ZIPCode) 

Monterey,  California  93943-5000 


8b  OFF'CE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBE 
(If  applicable)  I 


8c  ADDRESSfOty,  State,  and  ZIP  Code) 

1629  S.  Street 
Sacramento ,  CA  95S1- 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


vVOPt  UNIT 
ACCESSION  NO 


•3b  TIME  COVERED 
FROM  TO 


14  DATE  OF  REPORT  (Year,  Month.  Da))  '  S  PAL 

1989,  March  7  ] 


11  TITLE  (Include  Security  Classification ) 

COMPARISON  OF  MEASURED  AND  TRANSFORMED  DIRECTIONAL  WAVE  SPECTRA  USING 
LINEAR  .REFRACTION  MODEL  _  _  _  . _ 


12  PERSONAL  A'JTHORIS! 

Khalid,  Mo h arm  a  d  A . 


'  3a  type  Oc  REPORT 

Master's  Thesis 


16  SUPPLEMENTARY  NOtAt.'ON 

Tiie  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  offi¬ 
cial  reliev  or  position  of  tine  Denartment  of  Defense  or  the  U.S.  Government . _ _ 


18  SJBJFCT  TEPVS  {Continue  on  reverse  if  necessary  and  identify  by  b  cc1  number  » 

Ocean  Waves;  Directional  Spectra;  Wave 
Refraction 


19  A D S T  p A C 1  { Continue  on  reverse  if  necessary  and  identify  by  block  number)  __  n 

Deep  water  directional  wave  spectra,  measured  by  a  NDEC  3-meter  bucy 
off  Monte rev  Eav,  are  transformed  to  shallow  water  usinc  the  linear 
refraction  model  by  Dobson  (1967).  The  transformed  directional  spectra 
are  compared  with  measured  spectra  using  pressure  gauge  array  in  shallow 
water  at  Marina  and  Santa  Cruz.  Refraction  (Kr) ,  shoaling  (Ks)  and 
Jacobian  (J)  transfer  functions  are  computed. 

The  modeling  results  and  measured  wave  data  at  Marina  are  in  good 
agreement  for  eastcrlv  waves  at  hioh  frequencies.  On  the  other  hanc. ,  ^.he 
linear  refraction  model  gives  highly  variable  values  of  wave  solutions  at 
Santa  Cruz.  It  is  believed  that  the  unrealistic  wave  results  at  Santa 
Cruz  are  associated  with  the  complex  bathymetry  of  Monterey  Bay.  The 
present  work  examines  the  accuracy  and  limitations  of  linear  refraction 


ihRr-rvations  . 


:ol  hr*  F  i  old  ’ 


20  DSTPR  *iqn  A.A'iAB'li'Y  'if  ABSTRACT 

C3  jnc.  x;te d  unl'Mi-ed  □  save  as  »r> 


:ia  NAVE  Or  RESPONSIBLE  INDIVIDUAL 


21  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassif iea 


Thornton 


(408 


t:  ,lnrUidr  A  Rr)  , 

646-2847 


Code  6£Tm 


DD  Form  1473,  JUN  86 


Previous  editions  are  obsolete 

S/N  OU'2-LF-Ol  4-660  3 


_S E C i  ’ R I T_y  CLASJrUCA  "o’. 

UNCLASSIFIED 


l 


Approved  for  public  release;  distribution  is  unlimited 


Comparison  of  Measured  and  Transformed  Directional 
Wave  Spectra  Using  Linear  Refraction  Model 

by 

Mohammad  Asmatullah  Khalid 
Captain,  Pakistan  Navy 
B.Sc.  (Hon.)/  University  of  Karachi,  1981 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degrees  of 


MASTER  OF  SCIENCE  IN  HYDROGRAPHIC  SCIENCES 

and 

MASTER  OF  SCIENCE  IN  PHYSICAL  OCEANOGRAPHY 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1989 


Author ; 


Mohammad  A.  Khalid 


Approved  by: 


~V7 


Edward  B.  Thornton,  Thesis  Advisor 


Curtis  A.  Collins,  Chairman 
Department  of  Oceanography 


Gordon  E.  Schacher, 

Dean  of  Science  and  Engineering 


11 


ABSTRACT 


Deep  water  directional  wave  spectra,  measured  by  a  NDBC 
3-meter  buoy  off  Monterey  Bay,  are  transformed  to  shallow 
water  using  the  linear  refraction  model  by  Dobson  (1967). 
The  transformed  directional  spectra  are  compared  with 
measured  spectra  using  pressure  gauge  arrays  in  shallow 
water  at  Marina  and  Santa  Cruz.  Refraction  (Kr) ,  shoaling 
(Ks)  and  Jacobian  (J)  transfer  functions  are  computed. 

The  modeling  results  and  measured  wave  data  at  Marina 
are  in  good  agreement  for  easterly  waves  at  high  frequen¬ 
cies.  On  the  other  hand,  the  linear  refraction  model  gives 
highly  variable  values  of  wave  solutions  at  Santa  Cruz.  It 
is  believed  that  the  unrealistic  wave  results  at  Santa  Cruz 
are  associated  with  the  complex  bathymetry  of  Monterey  Bay. 
The  present  work  examines  the  accuracy  and  limitations  of 
linear  refraction  model  by  field  observations. 
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I. 


INTRODUCTION 


Wave  climatology  plays  an  important  role  in  coastal 
dynamics  as  most  energy  is  put  into  coastal  region  by  wind¬ 
generated  waves.  Initially,  ocean  wave  analysis  was  limited 
to  a  single  wave  train  defined  by  a  wave  height  and  a  period 
in  a  one-dimensional  spectrum,  but  in  the  last  couple  of 
decades  the  emphasis  has  been  placed  on  two-dimensional 
spectra . 

The  shoaling  waves  undergo  refraction,  reflection  and 
diffraction  mostly  influenced  by  the  wave  approach  angle  and 
the  bathymetry  of  the  area.  A  spectral  approac.i  requires 
the  superposition  of  a  large  number  of  progressive  sinusoi¬ 
dal  component  waves.  Therefore,  waves  of  small  amplitude 
(linear  wave)  are  assumed.  This  assumption  neglects 
reflection,  diffraction  and  postulates  long-crested  waves. 

The  area  of  Monterey  Bay  has  been  selected  for  this 
comparative  study  of  ocean  waves.  This  area  presents  an 
interesting  problem  for  predicting  directional  wave  spectra 
because  of  its  variable  shoreline  orientation  relative  to 
prevalent  eastward  waves  and  the  presence  of  the  Monterey 
submarine  canyon  bisecting  the  bay.  The  Monterey  submarine 
canyon  is  the  largest  in  the  western  hemisphere  and  causes 
substantial  perturbations  on  the  nearshore  wave  field. 


x 


Historically,  Monterey  Bay  was  one  of  the  first,  and 
presumably  the  most  intensively  studied,  locations  for  wave 
refraction.  Wave  refraction  diagrams  of  Monterey  Bay  were 
first  constructed  in  1948  by  Johnson,  O'Brien  and  Isaac 
(Wiegel,  1964).  The  first  intensive  refraction  studies  were 
performed  for  amphibious  landing  exercises  at  Fort  Ord 
(1954)  followed  by  studies  by  the  Corps  of  Engineers  (1958) 
for  small  craft  harbor  at  Santa  Cruz.  Wave  refraction  and 
littoral  drift  calculations  have  been  performed  by  Dorman 
(1969)  and  Arnal  et  al.  (1973).  All  these  studies  were 
performed  using  hand-drawn  refraction  diagrams.  More 
recently,  numerical  refraction  calculations  have  been  made 
for  a  number  oi  specific  marine  construction  sites  (e.g., 
Thornton  (1978,  1980))  and  erosion  studies  (Oradiwe,  1986; 
McGee ,  1986)  . 

The  objective  of  this  study  is  to  compare  ocean  wave 
measurements  with  the  linear  refraction  model.  The  deep 
water  two-dimensional  spectrum  is  transformed  to  shallow 
water  spectrum  using  transfer  function  computed  from  a 
linear  wave  refraction  model  by  Dcbson  (1967).  The  deep 
water  spectrum,  S0(fo,’g)  is  determined  from  data  acquired 
from  the  National  Data  Buoy  Center  (NDBC)  using  the  Fourier 
analysis  technique  as  proposed  by  Longuet-Higgins  (1963)  and 
Cartwright  (1963).  The  transformed  spectrum  S(f,  >),  is 
compared  with  measured  shallow  wave  spectra  for  which  the 
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data  had  been  acquired  from  pressure  gauge  arrays  located  at 
Marina  and  Santa  Cruz. 

The  background  theory  is  described  in  Chapter  II.  The 
data  acquisition  and  equipment  are  mentioned  in  Chapter  III. 
Chapters  IV  and  V  contain  discussions  and  the  conclusions. 


II.  THEORETICAL  BACKGROUND 


A.  WAVE  REFRACTION 

1 .  General 

The  refraction  of  ocean  surface  waves  on  the 
continental  shelf  is  necessary  for  forecasting  ocean  waves 
in  shallow  water.  Practical  techniques  for  calculating 
refraction  diagrams  were  developed  by  Johnson  (1948)  and 
Munk  and  Arthur  (1952).  The  graphical  technique  developed 
by  these  investigators  is  called  the  Orthogonal  Technique. 
With  the  advance  of  digital  computers,  the  numerical  methods 
were  introduced  by  Griswold  (1963),  Wilson  (1966)  and  Dobson 
(1967).  The  numerical  calculations  of  wave  orthogonals  (or 
rays)  can  be  quickly  and  accurately  plotted  for  engineering 
purposes.  The  study  of  wave  transformation  that  includes 
refraction  over  complex  underwater  topography  and  in  large 
areas  is  more  commonly  accomplished  through  numerical 
models.  (Hayes,  1980) 

2 .  Refraction  Theory 

In  this  study,  linear  wave  theory  is  applied  to 
shoaling  and  refracting  waves,  which  neglects  reflection  and 
diffraction  assuming  energy  propagates  only  in  the  direction 
in  which  the  waves  travel.  It  also  neglects  amplitude 
dispersion.  ine  wave  refraction  is  based  on  generalization 
of  Snell's  lav:  for  phase  speed  variation  and  conservation  of 
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energy.  The  phase  speed  of  the  wave  depends  upon  wave 
period  (or  frequency)  and  the  depth  of  water  at  particular 
locations.  Wave  energy  is  conserved  between  wave  orthog- 
onals,  where  orthogonal  is  perpendicular  to  the  wave  crest. 
The  orthogonal  indicates  the  direction  in  which  shoaling 
wave  progresses.  Based  on  linear  wave  theory,  it  is  assumed 
that  the  wave  period  will  remain  constant  but  the 
wavelength,  the  velocity  and  the  height,  will  vary  as  the 
wave  shoals  toward  the  beach. 

In  general,  the  phase  speed  of  a  small-amplitude 
linear  water  wave  is  expressed  as 

C2  =  ^  tanh(kd)  (2.1) 

where  g  is  acceleration  of  gravity,  k  the  wave  number  and  d 
is  the  depth  of  water  at  a  particular  location.  Consequent¬ 
ly,  in  shallow  water  the  phase  speed  of  the  wave  varies  with 
depth  only  and  Snell's  law  becomes  applicable  in  the  case  of 
straight  and  parallel  contours  as 
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where  is  an  angle  between  the  wave  crest  and  local  depth 
contours  and  C  is  phase  speed.  The  subscript  zero  refers  to 
deep  water  conditions. 
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The  energy  of  a  wave  when  it  reaches  the  beach  is  a 
function  of  the  wave  height,  and  the  wave  height  at  the 
shore  is  related  to  the  wave  height  in  deep  water.  The 
variation  of  wave  height  within  two  wave  rays  is  described 
(Mei,  1980)  as 


H  =  H 


0 


Cg01l/2  r^O.1/2 

Cg  J  1  b  J 


(2.3) 


or 

H  =  H0  KS  KR  (2.4) 

where 

Cc 

Ks  =  [_^]1/2  =  shoaling  Coefficient  (2.5) 

Kr  =  [-J1]-1-/2  =  Refraction  Coefficient  (2.6) 

where  b  is  the  distance  between  two  rays  and  Cg  is  group 
velocity  and  may  also  be  represented  in  the  form  of 

C  =  —  r i  +  2  kd  I 

3  2  sinh  2kd 

Using  the  relationship  C  =  C0  tanh(kd)  (Kinsman,  1965), 
Equation  (2.5)  may  be  written  as 
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(2 . 5a) 


k7 


k  -t  1/2 


1  + 


2  kd 


sinh  2kd 


-1/2 


From  Equation  (2.2),  the  change  in  angle  determines  the 
change  in  wave  crest  or  b  value.  Thus  the  value  of  KR  may 
be  computed  from 


b 

cos 


C^_ 

-‘0 


b 

COS  .1 


or 


Kr 


bg  =  rCOS  !  0,1/2 

*■  b  COS 


(2 . 6a) 


where 


.  -lr_ 
=  sin  iC 


sm 


C, 


In  any  linear  wave  refraction  model,  shoaling  and 
refraction  coefficients  are  first  computed  and  then  the  wave 
height  is  determined. 

3 .  Linear  Wave  Refraction  Model 

Munk  and  Arthur  (1952)  described  the  theory  of  wave 
intensity  along  a  refracted  wave.  Dobson  (1967)  followed 
the  theory  and  developed  a  linear  refraction  model.  The 
Dobson  model  is  briefly  described  here. 

The  definition  sketch  for  the  wave  ray  equation  is 
shown  in  Figure  2.1  in  that  is  tangential  to  circular 
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y 

A 


- ►  X 

Figure  2.1  Definition  Sketch  for  Wave  Ray 
Equations  (Dobson,  1967) 
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arcs  of  radius  C 't  and  (C+6C)  ft  and  S  and  n  are  very  small. 
Two  basic  characteristic  equations  for  the  propagation  of  a 
wave  ray  are  the  ray  curvature  equation 

Dj.  _  _  1  DC 

R  “  DS  C  Dn  (2.8) 

and  an  equation  for  the  ray  separation  factor  g  (Dobson, 
1967) 


d.-2  ,  d„-  DC  ,  D2c 

,2  +  dt  DS  n  2 

dt  Dn 


0 


(2.9) 


Dobson  manipulated  these  equations  using  Gauss 
integral  theory  into  a  convenient  form  where  Eqs.  (2.8)  and 
(2.9)  can  be  easily  applied  in  terms  of  cartesian 
coordinates  and  time  variable  as 


R  = 


,C 

—  -  cos 
X 


(2 . 10) 


p 


d_ 

dt 


+  a 


0 


(2.11) 


where 


P(t) 


2  [cos 


sin 


(2.12) 


q(t) 
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C  [  sin 


2  sin  cos 


x  v 


(2.13) 
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The  phase  speed  is  a  function  of  the  coordinates  so 
that  C  =  C(x,y).  The  wave  ray  can  be  constructed  with 
Equation  (2.10)  within  a  coordinate  system.  The  ray  curva¬ 
ture  R  is  computed  using  finite  difference  methods  between 
two  consecutive  points  along  a  wave  ray  by  iterating 

Ac.  =  R  A  S  (2.14) 


at  each  step.  So  the  new  coordinates  can  be  computed  as 

xs+  Q  =  XS  +  :S  COS  (2.15) 

Ys+iS  =  YS  +  1S  COS  7  (2*16) 


where : 


—  ^ 

K  =  -  ■  — 2-  ■ — - 

S  +  ‘‘S+AS 
2 

The  ray  path  is  constructed  grid  to  grid  over  the  sea  bottom 
until  it  reaches  the  required  location. 

The  ray  separation  factor  r  in  Equation  (2.11)  is 
computed  applying  Fox's  formula  for  a  second  order  linear 
equation.  This  is  a  finite  difference  method  and  requires 
that  the  separation  factor  be  known  at  the  starting  and 
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preceding  points.  This  condition  is  conveniently  satisfied 
by  starting  the  ray  in  the  deep  water  where  ,  is  equal  to  1. 

The  characteristics  of  the  wave  at  any  point  in  the 
coordinate  system  may  be  determined  from  the  solution  of 
Equations  (2.10)  and  (2.11)  provided  the  speed,  C(x,y)  is 
known.  However  it  is  common  to  consider  waves  of  different 
deep  water  characteristics  travelling  towards  the  coastline 
over  a  known  bathymetry,  so  that  it  is  more  convenient  to 
consider  the  problem  in  terms  of  water  depth,  d(x,y)  rather 
than  wave  speed.  Thus  these  equations  are  expressed  as 
functions  of  the  mean  water  depth. 

Dobson  used  a  quadratic  least  square  surface  to  fit 
depth  locally  at  intermediate  points  using 

d  =  i  +  ■  2x  +  •  3y  +  4x2  +  *'•  5xy  +  •  6y2  (2.17) 

where  the  partial  derivative  of  depth,  d,  with  respect  to  x 
and  y  is  expressed  as  functions  of  .  .  Therefore  the  local 
water  depth  can  be  obtained  by  interpolation  method. 

The  bathymetry  used  in  this  refraction  program  was 
obtained  from  the  original  NOAA  data  and  projected  onto  a 
six  second  (latitude  and  longitude  grid)  modified  Universal 
Transverse  Mercator  (UTM)  grid  after  screening  the  data  for 
errors  and  bad  points.  Insufficient  data  were  available  for 
off-shore  bathymetry.  Thus  it  required  interpolation  to 
fill  the  empty  grid  cells.  Piece-wise  linear  triangulation 
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was  used  to  interpolate  to  a  200  meter  rectangular  grid.  A 
nine-point  weighted  linear  average  was  applied  to  smooth 
bathymetry  (Oradiwe,  1986;  McGee,  1986). 

B.  DIRECTIONAL  WAVE  SPECTRUM 

1 .  General 

Directional  characteristics  of  ocean  waves  are 
important,  not  only  for  understanding  wave  generation 
mechanisms  but  also  for  practical  problems  such  as  wave 
forecasts  in  the  design  of  coast  and  harbor  protection 
works.  Directionality  is  equally  important  in  the 
refraction  of  ocean  waves.  Various  techniques  incorporating 
floating  buoys,  pressure  gauge  arrays,  electromagnetic 
current  meters  and  SAR  have  been  utilized  to  determine 
directional  properties  of  the  ocean  waves.  In  this  study 
data  have  been  acquired  from  the  NDBC  heave/pitch/roll  buoy 
for  deep  water  spectra  and  pressure  gauge  arrays  for  shallow 
water.  The  theoretical  background  of  these  measurements  is 
presented  in  the  following  sections. 

2 .  Deep  Water  Wave  Spectrum 

Longuet-Higgins  (1963)  suggested  a  method  to  compute 
directional  wave  spectrum  from  water  surface  elevation  and 
two  orthogonal  components  of  wave  surface  slope.  The 
surface  elevation  is  supposed  to  be  the  result  of  the 
superposition  of  a  large  number  of  long-crested  progressive 
sinusoidal  component  waves  in  random  phase  (Battjes,  1972). 
Several  other  investigators  (e.g.,  Cartwright,  1963;  Long, 
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1980;  Mitsuyasu,  1975;  and  Long  and  Hasselmann,  1979)  have 
suggested  and  refined  the  technique. 


The  directional  spectrum  S(f,6)  as  given  by  Longuet- 
Higgins  (1963)  is  expressed  in  the  form 

S(f,-0  =  S(f)  •  D  ( f , )  (2.18) 

where  S(f)  is  the  one-dimensional  frequency  spectrum  and 
D ( f ,  )  is  angular  distribution  function  with  the  condition 

2- 

D(f ,  >)d..  =  1  (2.19) 

0 

Equation  (2.19)  can  be  expanded  by  a  Fourier  series 

D  ( f ,  )  =  -  [i  -  ~  ancos(n)  +  bnsin(nO  ;]  (2.20) 

n=l 

The  coefficients  of  the  first  two  harmonics  can  be  estimated 
from  the  co-  and  quad-spectra  from  heave  (index  1),  and 
slope  signals  (indices  (1,2),  correspond  to  x,y  slopes) 

( Longuet-Higgins ,  1963;  Long,  1980)  as 

2 '  q12 

a  -j_  =  D  ( f  r  ■ )  cos  «  d  ,  (2.21) 


bi  = 


2 


sin  .  d 


(2.22) 


£] 


a2 


2- 


0 


D  (  f  ,  a  ) 


cos  2  a 


da 


C22  “  C33 


2 

K  C 
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(2.23) 


b2 


2- 

I"  D  (  f  ,  a  ) 

0' 


sin  2a  da 


2C 
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where  the  constant  K  is  defined  as 


(2.24) 


K  =  [C«  t 

L11 

A  heave/pitch/roll  buoy  and  pressure  sensor  array 
provide  a  measure  of  -  ,  r,x,  and  r.  y  at  fixed  points.  The  co- 
and  quad-spectra  of  pairs  between  these  three  quantities  are 
computed  using  the  Fast  Fourier  Transform. 

The  wave  directionalities  are  determined  by  the 
directional  distribution  function.  Equation  (2.20)  may  be 
rewritten  as 


D(f , 


JL 

2- 


r-^cos  (  - 


+  r2cos  2  ( a — 2 )  ] 


(2.25) 


where 


rx  =  [ a2  +  b2 ] 1/2  and  r2  =  [a2  +  b2]1/2 


In  the  case  of  unimodal  spectrum,  Longuet-Higgins 
(1963)  proposed  a  cosine  powered  distribution  function 
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D(f,a)  =  G  ( s) 


(2.26) 


cos 


^  (  u  -  a  ) 


i  2S 


where 


G(S)  = 


,2S-1 


(S  +  l) 
(2S+1) 


=  2 


2S-1  (S!) 


( 2  S  )  I 


is  a  normalizing  function  to  satisfy  Equation  (2.19).  The 
parameter  S  is  a  function  of  wave  frequency  and  a  is  the 
mean  direction  computed  from 

.  _  -1,Q13. 

=  tan  ( — )  =  tan  (- — ) 

al  ‘*-12 

Cartwright  (1963)  showed  that  the  value  of  S  is  determined 
by 


1  +  3r2  +  (r2  +  14r2  +  1)  1/2 
S2  =  2(1  -  r 2 ) 

The  cosine  distribution  model  is  used  most  in  uni- 
modal  wave  analysis.  The  same  expression  is  used  in  this 
study . 

3 .  Predicted  Shallow  Water  Wave  Spectrum 

The  shallow  water  wave  spectrum  can  be  predicted 
from  the  deep  water  spectrum  by  multiplying  the  squares  of 
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both  refraction  (Kp)  and  shoaling  (Ks)  coefficients  and  the 
Jacobian  of  the  direction  function  by  each  frequency 
component.  Hence  the  deep  water  spectrum  is  transformed 
into  shallow  water  spectrum  by 


2  2 

S(f  ,  )  =  S  o  ( f  o  ,  o )  Kr  Ks  J 


(2.27) 


where  the  Jacobian 


f(W 

dif,  .) 


Since  the  wave  period  remains  invariant,  J  only  depends  upon 
the  changes  of  local  wave  direction. 

4 .  Measured  Shallow  Water  Wave  Spectrum 

The  two-dimensional  shallow  water  wave  spectrum  is 
determined  from  the  data  acquired  from  shallow  water 
sensors.  The  same  technique  as  described  previously  for 
deep  water  spectrum  is  used.  In  this  study,  the  pressure 
gauge  arrays  have  been  utilized  in  which  subsurface  pressure 
fluctuations  are  measured  to  infer  the  surface  elevation 
(Kim,  1974;  Pawka  et  al . ,  1976). 

From  linear  wave  theory,  the  relation  between  wave 
elevation  (• )  and  subsurface  pressure  (P)  is  described  by 


ft)  =  i - H°vh  -k!  p(t) 

cosh  k ( a  +  z i 


(2.28) 
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where 


d  =  water  depth  in  meters, 

z  =  vertical  distance  from  mean  sea  level  to 
pressure  gauge  (z  is  positive  upward  from 
sea  level) , 

P  =  pressure  measured  in  meters  of  water  head 
about  a  zero-mean  pressure. 

The  term  in  the  brackets  in  Equation  (2.28)  is  the  spectral 
transfer  function  relating  pressure  to  surface  elevation  as 
a  function  of  frequency  and  depth.  Equation  (2.28)  can  be 
modified  to  include  differential  distance  and  elevation 
between  two  sensors  of  an  array  to  compute  co-  and  quad- 
spectra  of  surface  slope  as 


clx 


1  .  cosh  kd  1 2 

2  ^oosh  k(d+z)  J  IP 

A 


1  ,  cosh  kc  .  c  (4r) 

x  lcosh  k(a+z)  T,  ,1P 


(2.29) 


(2.30) 


where  P>;  is  subsurface  pressure  difference  between  two 
sensors  on  the  >:-axis  and  x  is  differential  distance  along 
the  same  axis.  Similarly,  pressure  power  spectra  and  co  and 
quad-spectra  can  be  computed  along  y-axis.  Once  the  powrer 
spectra,  co,  and  quad-spectra  are  obtained,  the  directional 
spectral  density  may  be  determined  using  the  same  technique 
outlined  earlier  for  deep  water  directional  spectrum. 
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III. 


DATA  ACQUISITION  AND  SENSORS 


A.  LOCATION  AND  TIME 

The  area  of  Monterey  Bay  was  selected  as  the  study  area 
which  poses  complex  problems  due  to  the  highly  irregular 
bathymetry  of  the  Monterey  submarine  canyon.  Three  sensors 
were  utilized  for  acquisition  of  wave  data.  The  locations 
of  the  sensors  are  as  follows  (see  Figure  3.1  for  bathymetry 
and  location) . 


No. 

Tvoe  of  Sensor 

Location 

Depth  CM) 

1 . 

NDBC  roll/pitch 
wave  buoy 

36°48 ' 00 ' 'N 
122°23 ' 59 ' 'W 

2003 

2  . 

Array  of  pressure 
gauges  (4) 

Marina 

3  6°4  2 . O'N 

12 1°4  8 . 9  '  W 

15.0 

3  . 

Array  of  pressure 
gauges  (4) 

Santa  Cruz 

3  6°57 .O'N 

1 2  2 °0 0 . 2  'W 

13 . 1 

Selected  wave  cases  from  January  1988  are  used  in  the 
analysis.  During  the  winter  season,  weather  conditions  are 
generally  favorable  for  creating  higher  energy  waves.  In 
addition,  an  abnormal  storm  was  encountered  in  the  region. 
To  select  a  particular  date  and  time,  the  conditions  for 
narrow  band  wave  spectrum  were  emphasized  to  simplify  the 
interpretations . 
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In  comparing  coincident  data  sets,  on  the  average  a  time 
lag  of  one  hour  was  assumed  for  the  waves  to  travel  from  the 
deep  water  buoy  to  the  shallow  water  array  sites. 

B.  DEEP  WATER  WAVE  DATA 

Deep  water  wave  data  for  the  pitch/roll  buoy  were 
obtained  from  the  NDBC.  The  NDBC  3-meter  buoy  measures 
heave,  pitch  and  roll  movements  with  respect  to  3-axes 
magnetometer.  The  magnetometer  axes  are  aligned  with  fore 
and  after  (bow  positive) ,  athwartship  (starboard  positive) 
and  vertical  (positive  upward)  axes  of  a  buoy.  The  analog 
sensor  measures  vertical  acceleration  and  displacement  in 
addition  to  hull  pitch  and  roll.  The  magnetometer  bow  and 
starboard  output,  together  with  vertical  displacement,  pitch 
and  roll  output  from  analog  sensors  are  sent  to  the 
Directional  Wave  Data  Analyzer  ( DWDA )  after  voltage 
amplification  and  filtering. 

The  hull  azimuth  angle  is  calculated  from  bow  and  star¬ 
board  magnetic  components  with  necessary  corrections  made 
for  pitch,  roll  and  hull  magnetic  fields.  The  magnetic  azi¬ 
muth  is  converted  into  true  azimuth  using  the  magnetic  vari¬ 
ation  for  the  site.  The  true  azimuth  is  then  used  to 
calculate  East-West  and  North-South  components  of  hull 
slope.  One  hertz  time  record  of  hull  displacement  and 
slopes  are  stored  in  DWDA  for  spectral  analysis  and  trans¬ 
mission  to  shore  station  through  GOES  (Steele  et  al.,  1985). 
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Data  provided  on  magnetic  tape  included  environmental 
parameters,  wave  spectra,  directional  wave  parameters  and  co 
and  quad  spectra.  Directional  wave  parameter  data  were 
utilized  to  determine  directional  wave  spectrum.  The 
unimodality  of  the  directional  waves  was  confirmed  by 
studying  the  regional  weather  maps  to  conveniently  apply 
cosine  angle  distribution  function.  Two-dimensional  spectra 
are  shown  in  Figure  3.2. 

C.  SHALLOW  WATER  WAVE  DATA 

The  calibrated  data  for  arrays  of  pressure  gauges  were 
obtained  from  the  Scripps  Institute  of  Oceanography.  The 
array  consists  of  four  pressure  gauges  located  in  a  six- 
meter  square  area  which  measures  pressure  in  centimeters  of 
water.  Each  sensor  is  located  about  50  to  100  cm  above  the 
bottom.  Pressure  power  spectrum  is  calculated  which  in  turn 
is  transformed  into  surface  elevation  and  surface  slopes 
x  and  y,  spectra  by  applying  linear  wave  theory  transfer 
function  as  described  in  Chapter  II. 

For  the  computation  of  the  cross-spectra,  the  following 
values  were  utilized: 

-  sampling  interval  =  -t  =  2  sec, 

-  total  data  points  =  2048, 

-  data  points  for  one  sample  =  128, 

-  total  record  length  =  32  min, 

-  Nyquist  frequency  =  0.25  Hz 
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.00391  Hz, 


-  frequency  bandwidth  =  if  = 

-  degrees  of  freedom  =  v  =  32. 

The  spectra  were  cut  at  frequency  of  0.2  Hz  in  the 
higher  frequency  range  due  to  spurious  overestimation  of  the 
transfer  function.  The  wave  period  for  high  quality 
estimate  is  approximately  5-20  seconds. 

To  ensure  accuracy  and  to  check  computations,  the 
results  were  compared  and  verified  with  the  published  data 
in  the  Coastal  Data  Information  Program  of  January,  1988. 
Two-dimensional  spectrum  at  the  peak  energy  frequency  is 
depicted  in  Figure  3.3  for  the  Marina  array  on  29  January 
1988  . 

D.  OCEAN  WAVE  TRANSFORMATION  PARAMETERS 

The  Dobson  linear  refraction  model  has  been  used  to 
estimate  the  shoaling  transformation  of  the  wave  direction¬ 
al  spectrum.  The  reverse  projection  method  is  employed. 
That  is,  the  refraction  program  was  run  backward  from 
shallow  to  deep  water,  commencing  at  a  shallow  water  loca¬ 
tion  in  required  depth,  for  Marina  15  meters  and  for  Santa 
Cruz  13.1  meters.  Rays  of  a  particular  frequency  were 
propagated  offshore  in  increments  of  .2  degrees  (0.2°)  over 
the  range  of  all  possible  incoming  wave  angles.  Once  the 
depth  is  deep,  the  rays  are  stopped  and  deep  water  wave 
directions  are  recorded.  The  rays  were  then  returned  to 
retrace  the  same  path  to  verify  the  backward  refraction  and 
to  calculate  transfer  coefficients  KR,  Ks  and  J. 
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MEASURED  SHALLOW  WAVE  SPECTRUM-MARINA 


To  include  refraction  over  the  shelf  to  the  north  of 


Santa  Cruz,  supplemental  bathymetric  data  were  meshed  into 
the  model.  The  shelf  was  more  coarsely  gridded  to  a  1  Km 
grid  because  of  the  large  area  involved  and  due  to  the  fact 
that  the  waves  are  primarily  in  deeper  water  while 
traversing  the  shelf.  The  transfer  parameters  calculated 
from  the  two  areas  were  multiplied  together  to  determine  a 
close  relationship  for  use  in  spectral  transformation.  Some 
example  ray  traces  are  shown  in  Figure  3.4. 

1 •  Shoaling  Coefficient 

2 

Shoaling  coefficients,  Ks ,  computed  from  the  model 
for  Marina  and  Santa  Cruz  sites  are  depicted  in  Table  3.1. 


TABLE  3 . 1 

SHOALING  COEFFICIENTS,  Kg 


Frea. 

Shoaling 

Coefficients 

(HZ) 

Marina 

Santa  Cruz 

(15  m) 

(13.1  m) 

0.030 

2.204 

2 .351 

0.040 

1 . 689 

1.796 

0 . 050 

1.389 

1 . 472 

0 . 060 

1.197 

1.263 

0 . 070 

1 . 068 

1.121 

0 . 080 

0.979 

1.022 

0 . 090 

0.918 

0.951 

0. 100 

0 .877 

0.902 

0.110 

0.851 

0.868 

0.120 

0.837 

0.847 

0..  13  0 

0.834 

0.836 

0.140 

0.839 

0.834 

0 . 150 

0.850 

0.839 

0 . 160 

0.867 

0 .850 

0 . 170 

0.887 

0.865 
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Figure 
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X 

(Continued) 


2 .  Refraction  Coefficient 

The  model  computes  refraction  coefficient,  KR,  for 
every  0.2°  of  nearshore  angle.  KR  values  were  averaged  for 
every  degree  of  deep  water  direction.  In  the  linear  refrac¬ 
tion  model,  an  assumption  is  made  that  the  underwater  topog¬ 
raphy  changes  are  small  and  gradual  but,  occasionally,  this 
assumption  does  not  hold  in  Monterey  Bay  because  of  the 
canyon.  As  a  result,  the  refraction  model  sometimes  gives 
erratic  estimates  of  the  KR  values.  Therefore,  all  KR 
values  above  three  were  truncated  to  the  value  of  three.  It 
is  also  observed  that  refraction  estimates  are  more  stable 
in  higher  frequencies  (short  period  waves)  than  lower 
frequencies  (long  period  waves).  See  Figure  3.5. 

3 .  Jacobian  (J) 

d  l0 

The  Jacobian  (-g— )  is  calculated  from  the  change  in 
the  deep  water  direction  as  a  function  of  near  shore  angles 
for  each  wave  frequency  for  all  possible  wave  ray  direc¬ 
tions.  A  plot  of  q  versus  'S  (Figure  3.6)  shows  a  slope  of 
nearly  45°  over  most  of  the  angles  close  to  the  beach 
normal.  The  slope  tends  to  increase  as  the  wave  ray 
undergoes  greater  refraction.  The  Jacobian  calculated  from 
these  values  (Figure  3.6)  gives  similar  results,  i.e.,  the 
Jacobian  estimates  are  constant  and  closer  to  one  when  the 
approach  angle  is  within  25°-30°  of  the  beach  normal  and 
increasingly  unstable/erratic  at  extreme  limits.  Values  at 
the  limits  may  not  be  useful  for  computation  of  directional 
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REFRACTION  COEFFICIENT  -  MARINA 


Figure  3.5  Refraction  Coefficients  for  Waves  at  Marina 


Figure  3.6  Plot  of  Deep  Water  Approach  Directions  vs  Ashore 
Angles  for  Waves  Arriving  at  Marina  Location 


spectra  particularly  in  the  lower  frequency  range  of  .07  Hz 
and  less. 


The  Jacobian  appears  to  behave  similar  to  the 
refraction  coefficients.  The  Jacobian  estimates  at  higher 
frequencies  are  more  stable  than  at  lower  frequencies. 
Figure  3.7  shows  the  Jacobian  estimates  at  different  fre¬ 
quencies  for  the  Marina  array  location. 

4 .  Transformation  of  Deep  Wave  Spectrum 

The  deep  wave  spectrum  is  transformed  to  shallow 
water  wave  spectrum  using  Equation  (2.27).  Transformation 
parameters  calculated  from  the  refraction  model  are 
multiplied  with  deep  water  directional  spectrum  to  obtain 
shallow  wave  directional  spectrum.  Figure  3.8  depicts  the 
transformation  of  0.09  Hz  wave  for  different  wave 
directions.  The  graph  of  Kg  is  not  shown  since  at  a 
particular  frequency  it  has  a  constant  value  for  all 


direct! ons . 


Figure  3.7  Jacobian  (J)  at  Various  Frequencies  for  Marina 


Figure  3.8  Transformation  of  Spectrum  from  Deep  Water 
to  Shallow  Water  at  Marina  for  Peak 
Frequency  0.09  Hz 
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IV.  RESULTS  AND  DISCUSSION 


A.  REFRACTION  MODEL  PREDICTIONS 

Transformation  parameters,  KR  and  Jacobian,  determined 
by  the  linear  wave  refraction  model  were  variable,  with  a 
tendency  to  increased  variability  with  increased  refraction. 
The  lower  frequency  waves  are  more  refracted  than  high 
frequency  waves.  Therefore,  these  coefficients  tend  to  vary 
smoothly  in  higher  frequency  (short  wave  period)  and 
increased  in  variation  with  a  decrease  in  frequency.  The 
variability  of  KR  and  J  is  greater  for  the  Santa  Cruz  site 
for  which  incident  waves  must  undergo  significant  refrac¬ 
tion  (see  Figures  4.1  and  4.2).  The  characteristics  of  each 
location  and  refraction  model  predictions  are  described 
here . 

The  Marina  location  is  shadowed  by  Point  Pino  in  the 
Southwest  and  Point  Santa  Cruz  in  the  Northwest.  High  fre¬ 
quency  waves,  0.11  Hz  and  above,  are  limited  to  directions 
from  235°-310°  at  the  array  site,  whereas  low  frequency 
waves  may  arrive  from  even  southerly  directions  after 
considerable  refraction  (refer  to  Figure  3.4).  Thus,  the 
overall  cut-off  limits  are  235°  in  the  South  and  310°  in  the 
North.  The  wave  rays,  at  low  frequency,  are  sparse  and  tend 
to  form  caustics  which  result  in  higher  variability  for 
estimates  of  coefficients.  The  variation  of  J  also  depends 
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REFRRCTION  COEFFICIENT 

( FREQUENCY  -  0.11  HZ) 
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Figure  4.1  Comparison  of  Refraction  Coefficients; 
Marina  and  Santa  Cruz 
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Figure  4.2  Comparison  of  Jacobian  at  0.13  Hz 
at  Marina  and  Santa  Cruz 
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upon  the  direction  of  approaching  waves  in  addition  to  the 
frequency.  If  the  wave  ray  coincides  with,  or  is  closer  to 
the  beach  normal,  the  J  values  smoothly  vary  over  the 
directions.  As  the  angle  of  approach  increases,  J  tends  to 
vary  with  increasingly  high  estimates. 

The  Santa  Cruz  site  poses  more  complex  problems  than 
Marina  because  of  its  location.  The  waves  approaching  from 
the  westerly  direction  makes  almost  a  90°  turn  to  arrive  at 
the  shoreline.  The  low  frequency  waves  travelling  north  and 
north-eastward  may  arrive  at  this  location.  The  ray  pattern 
shows  (Figure  4.3  a,b)  scattered  and  crossing  rays  (caus¬ 
tics)  with  occasional  or  a  few  wave  rays  from  the  west  or 
higher  approach  angle  at  low  frequency.  This  gives  very 
large  and  highly  variable  estimates  of  J.  The  ray  pattern 
becomes  more  regular  with  increase  in  frequency.  It  is  also 
noticeable  from  the  ray  patterns  that  the  J  estimates  are 
generally  large  and  more  variable  compared  with  the  Marina 
values  (Figure  4.2).  The  results  of  KR  estimates  behave 
similarly . 

B.  WAVE  DATA  SELECTION 

Data  for  wave  measurements  for  the  month  of  January  1988 
were  examined.  A  number  of  cases  were  analyzed  out  of  which 
three  are  described  here.  The  selected  days  are  4  January 
0800  PST,  18  January  at  0200  PST  and  29  January  at  2000  PST. 
These  periods  provide  narrow  band  energy  density  at  differ¬ 
ent  frequencies  and  direction.  Figures  4.4,  4.5,  4.6  and 
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(a)  Low  Frequency 


Figure  4.3  Wave  Ray  Traces  at  Santa  Cruz 
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(b)  At  High  Frequency 
Figure  4.3  (CONTINUED) 
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MEASURED  SPECTRAL  DENSITY 
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Figure  4.5  Measured  Spectral  Density  for  29  January  1988, 

2000  PST 


MEASURED  SPECTRAL  DENSITY 


Tables  4.1  and  4.2  provide  a  comparison  of  measured 
directional  spectra  at  the  three  locations  for  each  day. 


TABLE  4.1 

DEEP  WATER  ENERGY  DENSITY  WITH  MEAN  DIRECTIONS 


FREQ 

3  JAN 

0800  PST 

18  JAN 

0200  PST 

29  JAN 

2000  PST 

(Hz) 

DIR 

ENERGY 

DIR 

ENERGY 

DIR 

ENERGY 

( 1/sec) 

(Deg) 

(m2-sec) 

(Deg) 

(m2-sec) 

(Deg) 

(m2-sec) 

0.0300 

173.4 

0.000 

283.4 

0.000 

183 . 3 

0.000 

0.0400 

163 . 8 

0.000 

297.3 

0.000 

234.3 

0.000 

0.0500 

240.5 

0.000 

265.0 

0.000 

223 . 7 

0.000 

0 . 0600 

214 . 2 

0.000 

271.5 

26.658 

224 . 1 

0.804 

0.0700 

256.6 

1.245 

284 . 6 

72.766 

256.5 

2.125 

0.0800 

270.0 

1.569 

292 . 8 

70.196 

273 . 5 

6.346 

0.0900 

268 . 5 

1.198 

302.2 

38 . 581 

276.6 

8 . 018 

0 . 1000 

242 . 8 

1.510 

311.8 

22 . 677 

277 . 5 

4.256 

0. 1100 

236.7 

2 . 398 

308.7 

12.638 

276.6 

3 . 007 

0.1200 

228 . 3 

2.755 

313.5 

11.789 

272.2 

2 . 329 

0.1300 

227.7 

2.881 

321.3 

10.475 

268 . 8 

1.570 

0.1400 

223.9 

2 . 190 

313 . 6 

8.089 

269.5 

1.113 

0.1500 

217.6 

1.699 

321.7 

3.894 

268.7 

0.675 

0.1600 

210.1 

1.418 

322.5 

2.687 

265.5 

0.609 

0.1700 

190.5 

1.347 

312.1 

2 .315 

259.3 

0.350 

0 . 1800 

179.2 

0.889 

309 . 6 

1.940 

259.6 

0.233 

0.1900 

175.8 

0.899 

317.2 

1.550 

267.2 

0.184 

0.2000 

173 . 6 

0.834 

311.2 

1.147 

245.1 

0.146 

0.2100 

165.2 

0.670 

313.8 

0.767 

235.6 

0.166 

TABLE  4 . 2 

WAVE  CHARACTERISTICS  OF  SELECTED  CASES 


CHARACTERISTICS 

3  JAN 

88 

18  JAN 

88 

2  9  JAN 

88 

^sign . 

2.2  m 

7.9  m 

2.3  m 

Freqpeak 

0.13 

Hz 

0.07 

Hz 

0.09 

Hz 

Frec2mean 

0.18 

Hz 

0.10 

Hz 

0.12 

Hz 

Directionmean 

230° 

300° 

275° 
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The  California  coast  was  hit  by  an  unusually  severe 
storm  during  16-18  January  1988  which  resulted  in  high 
waves,  strong  winds  and  sea  level  well  in  excess  of 
predicted  values  (Cayan,  et  al.,  1988).  Widespread  damage 
occurred  along  the  Southern  California  coast.  This  storm 
developed  about  300  NM  west  of  San  Francisco  on  16  January 
1988  and  traveled  south-east  at  about  33  knots.  The  storm 
center  passed  off  Monterey  Bay  on  the  morning  of  17  January 
1988  continuing  on  its  SE  track  and  making  landfall  at  Avila 
Beach  in  Central  California.  The  bay  experienced  record 
high  waves  14  hours  after  the  storm  passed  from  the  area.  A 
significant  wave  height  of  7.9  meters  was  recorded  on  the 
morning  of  18  January  1988.  The  dominant  wave  period  was 
about  14  sec  with  a  direction  of  284°  with  narrow  band 
energy  measured  at  deep  water  and  shallow  water  locations 
(Figure  4.4). 

Deep  water  energy  density  with  mean  direction  for 
selected  cases  are  given  in  Tables  4.1  and  4.2.  On  29 
January  1988  measured  spectral  density  of  deep  water  waves 
and  shallow  coastal  site  are  shown  in  Figure  4.5. 

Reviewing  Tables  4.1  and  4.2  and  other  published  data, 
deep  water  waves  generally  travel  eastward  during  the  month 
of  January  with  slight  variability  of  about  20°  on  either 
side  of  true  East.  This  direction  is  dominant  in  the 
energetic  middle  frequency  range  of  0.08  to  0.13  Hz.  At 
other  frequencies,  the  direction  tends  to  be  variable. 
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The  third  case  of  3  January  1988  provides  a  different 
view  from  the  normal  as  it  contains  sea  wave  energy  at  0.13 
Hz  frequency  with  a  mean  southerly  direction  of  228°  (see 
Figure  4.6).  In  this  case,  deep  water  waves  undergo  less 
refraction  to  arrive  at  the  Santa  Cruz  array. 

C.  COMPARISON  OF  MEASURED  AND  CALCULATED  SPECTRA 

The  directional  spectra  from  wave  measurement  data  were 
determined  using  a  cosine-powered  distribution  function 
(Equation  (2.26)).  The  value  of  the  S  parameter,  in  most 
cases,  was  less  than  20  and  relative  directional  energy  was 
in  a  broad  band  of  50°  to  60°  spread. 

Predicted  and  measured  shallow  wave  directional  spectra 
for  Marina  are  compared  in  Figures  4. 7-4. 9.  The  deep  water 
wave  approaching  directions  are  limited  between  235°  to 
310°.  These  limits  change  slightly  with  decrease  in  wave 
frequency.  In  the  low  frequency  range,  the  predicted 
spectral  energy  is  overestimated  (see  Figures  4.7  and  4.8). 
The  deep  water  spectra  of  high  frequency  waves  from  the  West 
can  be  transformed  into  shallow  water  with  reasonaoly  good 
estimation  comparable  with  actual  wave  measurements.  The 
predicted  spectrum  from  westerly  waves  is  shown  in  Figure 
4.9.  The  predicted  mean  shallow  water  direction  tended  to 
be  more  aligned  with  the  deep  water  spectrum  than  the 
measured  shallow  water  spectrum. 

The  transformation  does  not  appear  to  appropriately 
shift  the  mean  direction  of  spectral  peak  energy  of  deep 
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Figure  4.7  Comparison  of  Shallow  Water  Spectra  at  Marina 
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Figure  4.9  Comparison  of  Predicted  and  Observed  Shallow  Wave  Spectra 
at  Marina  for  a  Frequency  of  0.09  Hz  for  29  January  1988 


water 

wave  spectra.  Two 

additional 

cases  of 

5  January 

1988 

and  30  January  1988  are 

shown  in  Figures  4.10  and  4.11 

for 

high 

frequency  spectra. 

The  mean 

approach 

direction 

for 

these 

frequencies  in  deep 

water  and 

shallow  water  are: 

Deep 

Shallow 

5  January  1988 

0.12  Hz 

255° 

272° 

30  January  1988 

0.13  Hz 

292° 

281° 

The  predicted  spectrum  in  shallow  water  also  indicates  high 
spectral  density  in  the  mean  direction. 

At  higher  frequency  sea  waves,  the  refraction  model 
overpredicts  spectral  density  by  less  than  18%  when  summed 
over  all  angles  at  a  particular  frequency.  For  swell  waves, 
the  results  are  largely  exaggerated. 

Results  of  the  Santa  Cruz  location  are  depicted  in 
Figures  4.12  and  4.13.  The  results  are  quite  different  from 
Karina.  The  predicted  spectra  are  highly  overestimated  at 
lower  as  well  as  higher  frequencies.  The  relatively  simple 
case  of  3  January  1988  was  particularly  selected  for  Santa 
Cruz  to  consider  ocean  waves  from  the  southwest  direction. 
Even  this  simple  case  does  not  provide  satisfactory  results. 

These  differences  are  broadly  associated  with  bathy¬ 
metry.  The  non-linear  effects  due  to  diffraction  and 
currents  were  neglected  which  may  introduce  errors.  It  is 
surmised  that  the  complexity  of  bathymetry  in  Monterey  Bay 
may  contribute  significantly  to  overestimation  and  variation 
of  transformation  parameters. 
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sdicted  and  Observed  Shallow  Wave  Spectra 
requency  of  0.12  Hz,  5  January  1988 


SHALLOW.  WAVE  SPECTRUM  (0.13  HZ) 
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Figure  4.11  Comparison  of  Shallow  Wave  Spectra  at  Marina 
for  30  January  1988 
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Figure  4.12  Comparison  of  Shallow  Wave  Spectra  at 
Santa  Cruz  for  29  January  1988 
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V.  CONCLUSION 


Deep  water  directional  wave  spectrum  was  determined  from 
data  acquired  from  the  NDBC  buoy  and  transformed  to  shallow 
water  using  the  linear  wave  refraction  model.  The 
transformed  spectrum  was  compared  with  measured  shallow  wave 
spectra  at  Marina  and  Santa  Cruz. 

The  transfer  parameters  of  refraction  coefficients  and 
Jacobian  predicted  by  the  linear  refraction  model  are  highly 
dependent  upon  wave  period  and  direction  of  approaching 
waves.  The  Jacobian  appears  to  be  more  sensitive  to  wave 
period  and  direction  than  the  refraction  coefficient.  These 
coefficients  are  generally  stable  at  higher  frequencies 
and/or  when  the  wave  ray  is  normal  or  nearly  normal  to  the 
coastline . 

The  model  tends  to  overestimate  the  transformed 
spectrum.  The  estimates  are  less  than  18%  in  excess  of  the 
measured  spectral  density  for  easterly  waves  at  frequencies 
of  0.09  Hz  and  higher  for  the  Marina  site.  At  lower 
frequencies  the  results  are  largely  exaggerated. 

The  model  provides  unrealistic  large  variations  in  wave 
height  calculations  for  waves  at  Santa  Cruz.  This  is 
presumably  due  to  its  location,  orientation  of  shoreline 
with  prevailing  north-easterly  wave  direction  and  the 
presence  of  Monterey  submarine  canyon. 
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It  is  clear  that  the  strong  irregularities  of  bathymetry 
in  Monterey  Bay  causes  non-linear  effects.  The  linear 
refraction  model  is  not  suitable  for  handling  such  problems. 
Further  modeling  efforts  including  non-linear  effects  are 
needed  for  Monterey  Bay  wave  study. 
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